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Abstract: We review a formulation of a renormalization-group scheme for Hamiltonian systems with 
two degrees of freedom. We discuss the renormalization flow on the basis of the continued fraction 
expansion of the frequency. The goal of this approach is to understand universal scaling behavior of 
critical invariant tori. 
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1 Introduction 

Recently, renormalization-group ideas have been proposed to describe the breakup of invariant tori 
for Hamiltonian systems with two degrees of freedom ||, ||, ||, |l| . Most of the numerical work has 
been done for the golden mean torus (mainly for practical convenience). In this article, we extend 
the renormalization-group transformation to more general frequencies. 

We consider the following class of Hamiltonians with two degrees of freedom, quadratic in the actions 
A = (A 1 ,A 2 ): 

H(A, cp) = im(^)(0 • A) 2 + K + 9 (^)0] • A + /( £ ), 

where m, g, and / are even scalar functions of the angles ip = (ipi, Lp-i) € T 2 = [0, 27r) 2 ; the frequency 
vector of the considered torus lo — (oj, —1) satisfies a Diophantine condition \uj_ q ■ y\ > <j\v\~ t for 
all v_ G Z 2 , and for some a > and r > 1. The rank of the matrix d 2 H/dA 2 is equal to one (this 
matrix is proportional to the projection operator on the f2-direction) . Thus, det [d 2 H/dA?] = 0, 
i.e., the Hamiltonians we consider do not satisfy the twist condition. There is however a twist in a 
particular direction of the actions characterized by the vector VL = (l,ct), which is what is needed 
to obtain KAM stability of Diophantine tori for small perturbations . 
If one takes g and / equal to zero, the equations of motion are 

# = -5t<2-^. f =.»< £ )<2-4)2 + ^ 

Then ^4 = defines an invariant torus, and the motion on this torus is quasiperiodic with frequency 
vector uj_ q . For g and / sufficiently small, a KAM theorem Q proves the persistence of this torus for 
arbitrary m, provided that its mean value is nonzero. It is not necessary to have m small in order 
to prove the existence of the torus with frequency vector cj . This remark allows us to perform 
canonical transformations that stay within the class of Hamiltonians quadratic in the actions. 
As one increases the perturbation (consisting of g and /), the torus gets deformed. This deformation 
is related to families of nearby periodic orbits that accumulate at the torus. This accumulation 
motivates the setup of a renormalization transformation combining a rescaling of phase space and 
an elimination of the irrelevant part at each scale. 

An attractive (trivial) fixed point of the renormalization represents the phase where the torus exists. 
An hyperbolic (nontrivial) fixed point (or more generally an hyperbolic fixed set) corresponds to a 
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transition where the torus breaks up. 

For the golden mean torus with frequency vector lj = (7~\ —1) where 7 = (1 + v5)/2, numerical 
studies suggest that the critical surface is the stable manifold (of codimension 1) of a nontrivial fixed 
point. Also for quadratic irrationals (those with periodic continued fraction expansion), one expects 
a nontrivial fixed point for a certain renormalization operator. For nonperiodic continued fraction 
expansion, hyperbolic fixed sets instead of fixed points are conjectured [ll], [l^, [l4|, [lS|, p"5| |. 

2 Renormalization transformation 

We describe the renormalization scheme for a torus with arbitrary frequency vector cj = (u, — 1) 
where u> G]0, 1[. This renormalization relies upon the continued fraction expansion of to: 

iv = ^ = [a ,ax, .. .}. 

dQ H ; 

di H 

The best rational approximants are given by the truncations of this expansion p k / qk — [a-o , &i , ■ ■ ■ ,a k = 
00]. The corresponding periodic orbits with frequency vectors y_ k = {qu,Pk) (called "resonances" 
in what follows) accumulate at the invariant torus. This family of frequency vectors satisfies 
l^o ' ^fe+il < l^o ' ^fcl anc ^ li m feoo |w • v_ k \ — Q, as it can be seen from the relation 

v k = N ao ■ ■ ■ N ak _ 1 v , 

where = (1,0) and N ai denotes the matrix 

"«-(? 0)' 

Moreover, lu q ■ y_ k and tj Q ■ v_ k+1 are of opposite sign (as the stable eigenvalue of N ai is negative); 
thus the torus is approached from above and from below by the sequence of periodic orbits with 
frequency vectors {v k }- 

We notice that if the continued fraction expansion is periodic with period s, then one can extract 
from the sequence {v k }, s families of periodic orbits which accumulate geometrically, with the same 
ratio, at the torus. 

The word resonance refers to the fact that the small denominators tv ■ v_ k that appear in the 
perturbation expansion are the smallest ones, i.e., |w -u k \ < \lv ■ y\ for any v_ = (q,p) different from 
zero and v_ k , and such that \q\ < qu+i- 

The main scale of the perturbation is defined by Vq for the torus with frequency uj. We will denote 
this scale by [j£ ,ct;]. The next smaller scale is [i/j^o;]. The renormalization transformation changes 
the coordinates such that the next smaller scale becomes the main one, i.e., the main scale is now 
[z^ ,a/] where u>' = [0,1,0,2, ■ ■ ■]■ As the frequency is changed (the continued fraction expansion is 
shifted to the left), the sequence of resonances is mapped into the sequence 

j/ k =N ai ■■■N ak _ 1 v . 

The renormalization transformation is a map (m,g, f,uj,a) <— > (m! , g' , f ,lu' , a'), where a denotes 

the second component of f2 = (1, a). 

The transformation lZ aQ consists of four steps: 

(1) A shift of the resonances constructed by the condition y_ x 1— > v^: we require that cos[(ao, 1) • (f\ — 
cos[(l,0) • ip']. This change is done via a linear canonical transformation (A,ip) 1— > (A^A, N ao ip). 
This step changes the frequency vector uj into u/ Q — (u)', —1), since N aa uj Q = —u>u/ . 

(2) We rescale the energy (or time) by a factor uj~ x , and we change the sign of both phase space 
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coordinates (A, ip) > {—A, —ip), in order to have u/ as the new frequency vector, i.e., the average 
of the term linear in the actions is of the form u/ ■ A. Furthermore, O = (1, a) is changed into 
tf = {l,(a +a)- 1 ). 

(3) Then we perform a rescaling of the actions: H is changed into H(A,ip) = \H(A/\,tp) with 
A = X(H) such that the mean value of m is equal to 1, i.e., A = W (ao + a) 2 {m). 

(4) A canonical transformation that eliminates the nonresonant part of g and /. 

The choice of which part of the perturbation is resonant or not is somewhat arbitrary. Recall that 
what is relevant is the accumulation of resonances; therefore a possible choice of irrelevant modes 
could include all modes except the resonances. However, it is desirable from a numerical point of 
view not to eliminate too many modes. A reasonable choice is the set C of integer vectors v_ such 
that \v% | > | v\ | . We notice that the relation defining v_ k = (qk,Pk) shows that qk > Pk for fc > 0, 
i.e., the resonances are not elements of C. 

From the form of the eigenvectors of N a<J , one can see that every i/£Z 2 \ (0,0) goes into C after 

sufficiently many iterations of matrices N ai (as the unstable eigenvector of N^ 1 , which is w , points 

into C). In other terms, a resonant mode at some scale turns out to be a nonresonant one at a 

sufficiently smaller scale. We notice that (0,0) is not an element of C, i.e., it is resonant. 

We eliminate completely all the nonresonant modes of g and / by a canonical transformation, 

connected to the identity, which is defined by iterating KAM-type transformations. 

We denote I - the projection operator on the nonresonant part: 

and d the derivative with respect to the angles ip: d = d/dp. 

The KAM iterations we perform are generated by functions linear in the actions, and it allows us, 

following Thirring ]l7[ , to remain quadratic in the actions at each step. 

This can be seen by working with Lie transformations Us '■ (A, cp) i— » (A , ip') generated by 

S(A, ip) = Y(p)£l ■ A + Z(p) +aQ-ip, 

characterized by two scalar functions Y and Z, and a constant a. The expression of the Hamiltonian 
in the new variables is obtained by 

H' = Ho Us = e § H 

= H + {S,H} + {S,{S,H}}/2\ + --- , 

where { , } is the Poisson bracket of two functions of the action and angle variables 

* ' 1 ~ dip ' dA ~ OA ' dip ' 

From this equation, one can see that H' is again quadratic in the actions: OS /dp and dH/dA are 
linear in the actions, dS/dA is action-independent and dH/dip is quadratic; therefore {S, H} is 
quadratic in the actions, and in fact, quadratic in the f2 • A variable. 

The generating function is determined such that it eliminates the nonresonant modes of g and /. 
This cannot be defined directly, so one iterates an infinite number of steps of such transformations 
such that one iteration eliminates the order e while producing an order e 2 . Following Ref. [Q, we 
have the equations 

cj ■ dZ + I - / = const, 

ui -dY + l~g + I- (mn ■ dZ) + afi 2 I"m = 0, 
(mjafl 2 + (g) + (mQ • dZ) = 0. 
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The constant a corresponds to a translation in the action variables which has the purpose of elimi- 
nating the mean value of the linear term in the variable f2 • A. We notice that Y and Z only contain 
nonresonant modes, e.g., Y — I~Y, and that Z (resp. Y) is chosen to reduce I - / (resp. l~g) from 
e to e 2 . These equations are solved by Fourier series: 

Y(u>) = V \g v + (roQ ■ dZ) v + m„afl 2 } e^. 

— UJr. ■ V L - - - J 

v£C — - 

For each iteration, we express the Hamiltonian in the new action and angle variables. This is done 
recursively, following Ref. ||. This KAM-type iteration adds terms of order e to m and to the 
resonant part of / and g. 

We iterate this procedure in order to reduce completely the nonresonant part of / and g: 
H' = H o Uh, where Uh — C^Si ° Us 2 ° • • ■ Us„ ° • • ■ , 

where I~f = I~g = 0. We notice that step (4) does not change w and f2 [as opposed to steps (1), 
(2), and (3)]. 

The transformation Uh is rigorously defined for a sufficiently small perturbation, consisting of / and 
g (see Ref. but the convergence in the whole domain of existence of the torus is a conjecture 

based on numerical observations. 

In summary the renormalization transformation acts as follows: First, some of the resonant modes 
are turned into nonresonant ones (by a rescaling of phase space). Then a KAM-type iteration 
eliminates these nonresonant modes, while slightly changing the resonant ones. 

3 Comments 

The numerical implementation of the renormalization scheme for a given frequency to shows that 
there are two main domains separated by a surface: one where the iteration converges to / = g = 
and the other where it diverges to infinity. 

The conjecture is that the boundary of the domain of convergence of the transformation, dD, coin- 
cides (at least locally, not too far away from the nontrivial fixed set) with the critical surface (where 
the torus is critical). This is by no means trivial since the transformation is based on properties 
that are valid for small / and g. However, for a one-parameter family, the numerical evidence of 
the coincidence between the critical coupling (where the torus breaks up, determined by Greene's 
criterion) and the value of the parameter where the iteration starts to diverge, indicates that we can 
expect dD to coincide with the critical surface, at least in a large region containing the nontrivial 
fixed set. 

In this section, we analyze the renormalization flow on the basis of the continued fraction expansion 
of the frequency. 

Quadratic irrational frequencies- We start by analyzing the effect on u> and a of s renormalization 
steps. Denote by {bj} the continued fraction expansion of a: a = [bo, b\, . . . ]. The renormalization 
^•a s _i^a s _ 2 ' ' -T^aa changes u = [ao, ffli, ■ • •] into [a s ,a s +i, •■•], andainto [a s _i,a s _ 2 , • ■ ■ , a ,&o,&i, ■••]• 
If lo has a periodic continued fraction expansion of period s, i.e., to = [(ai, . . . , a s )oo], one expect 
to have a nontrivial fixed point on the critical surface of the renormalization transformation in 
which one step is defined by the composition 72- 0s 72- 0s _ 1 • ■ - TZ ai - We notice that a converges to 
[(a s , . . . , ai)^]. Therefore f2 converges to the unstable eigenvector of the matrix N as ■ ■ ■ N ai (the 
stable eigenvector of this matrix is co ). 

The nontrivial fixed point associated with lo defines a universality class that we characterize with 
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critical exponents such as the total rescaling of phase space (product of the s rescalings), and the 
unstable eigenvalue of the linearized map around the fixed point. 

Two frequencies u>\ and L02 having the same periodic tail (and different first entries) in their contin- 
ued fraction expansions, belong to the same universality class. The initial integers in the continued 
fraction expansion are irrelevant. 

Associated with this nontrivial fixed point, we also have nontrivial fixed sets related to the nontrivial 
fixed point by symmetries. Therefore these hyperbolic sets belong to the same universality class. 
According to Ref. H) > these sets are given by the intertwining relation 

7t ao 0%= Tff ao e o 1Z ao , 

where 7g is defined as the translation (Tgf)((p) — f(ip + 9). Applying the relation 

K as o • ■ • o TZ ai oTe_ = TN as -N ai e • • • K ai 
to the fixed point H*(tp), we have 

n a3 o---oK ai H*{tp + 9) = H*(ip + N aa ---N ai 9). 

The map 

9^N a3 ■■■N ai e mod 2tt, 
gives the nature of the orbits to which the transformation converges. 

For instance, for the golden mean case, one has a nontrivial fixed point and a nontrivial 3-cycle in 
the space of even Hamiltonians. For Hamiltonians without parity restriction, the nontrivial fixed 

sets can be labeled by the orbits of the Anosov map TV = f J J ~\ (N 2 is Arnold's cat map) . For 

u> = [(1, 2)oo], one has a nontrivial 2-cycle and two fixed points in the space of even Hamiltonians. 
In the perturbative regime, there exists a geometrical accumulation of a certain sequence of periodic 
orbits with a "trivial" ratio (trivial in the sense that it is explicit, see comment in the previous 
section). The fact that it also happens in the critical case, but with a nontrivial ratio (which is 
the unstable eigenvalue of the linearized map around the nontrivial fixed point), implies universal 
self-similar properties of the critical torus |)| [l6] . 

As already mentioned in Ref. EJ, lZ ai becomes more and more singular as 04 increases. Therefore 
the picture we present is valid only for bounded ai. For large m, there is, to our knowledge, no 
description of the renormalization flow and its properties, even from a numerical point of view. 

Nonquadratic irrational frequencies- In that case, we cannot expect any geometrical accumu- 
lation of periodic orbits even if the Hamiltonian is close to integrable (perturbative regime). The 
renormalization has no fixed point, but instead its flow is related to a chaotic trajectory of the Gauss 
map 

LO I > 1/u — [1/lo], 

where [] denotes the integer part. We notice that a follows the inverse of this trajectory. 
For a typical w S]0, 1[, it is numerically impossible to figure out what could be the critical set, as 
one only works with finite precision (and thus a finite number of entries in the continued fraction, 
determined by iterating the Gauss map). 

Instead, one should give a sufficient number of entries in the continued fraction expansion, and then 
iterate the renormalization transformation on the critical surface. These entries should be bounded 
in order to avoid singularities mentioned above. 

On the critical surface, one can conjecture the existence of a critical strange attractor. But these 
hyperbolic sets should be conceptually different from the sets found in the quadratic irrational case 
(related to the nontrivial fixed point by symmetries): there is a continuous distribution of critical 
exponents (e.g., rescalings). The mean-rescaling (geometric mean value) and the largest Lyapunov 
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exponent characterize the universal class associated with u>. 

In order to have information about a typical lu g]0, 1[, an ergodic renormalization for random con- 
tinued fractions has been proposed jl^, ^5|. It consists in determining universal parameters by 
averaging over a large number of random continued fractions, constructed following a given proba- 
bility distribution for the coefficients. 

The conjecture is that ergodic renormalization trajectories converge to a strange attractor (which 
contains all the nontrivial fixed sets obtained for quadratic irrationals). Lyapunov exponents and 
other quantities such as the mean-rescaling are universal. Related ideas have been proposed for 
circle maps in Refs. [l2f . 

Extension to higher dimensional systems— We have seen that the existence of a nontrivial fixed 
point was based on a geometrical accumulation of a certain sequence of periodic orbits. These were 
given by the truncations of the continued fraction expansion (for quadratic irrationals). For three 
degrees of freedom, we lack of a theory that generalizes the continued fractions. However, we can 
choose a torus such that it is a geometrical accumulation of periodics orbits, i.e., the sequence of 
resonances is generated by a single matrix (as it was the case for the golden mean). It allows us 
to define a renormalization-group transformation with a fixed frequency vector. Even in that case, 
preliminary studies || suggest that one can expect a strange attractor instead of a nontrivial fixed 
point. This feature depends strongly on the spectrum of the matrix which generates the resonances. 

4 Analyticity properties of the nontrivial fixed point 

As the renormalization transformation described in the previous sections does not reduce the non- 
resonant part of m (it only reduces g and /) in order to remain quadratic in the actions, the shift 
of the Fourier modes [step (1)] reduces the analyticity of to at each step. As a consequence, the 
nontrivial fixed point is expected to have a nonanalytic to. The numerical implementation of the 
transformation gives accurate computation of the critical exponents characterizing the universality 
class, despite the fact that to is nonanalytic at the nontrivial fixed point. At present, this fact is 
not well understood. 

A method that can yield an analytic nontrivial fixed point is to eliminate the nonresonant part of 
to together with the one of g and / Q. The main drawback of this method is that one cannot 
remain quadratic in the actions. This drastically complicates the numerical implementation of the 
transformation as one needs to work with N scalar functions instead of 3, N denoting the numerical 
truncation of the Taylor series in the actions. 
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